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Abstract: This paper briefly reports some results of our study on the application 
of a decentralized adaptive control approach to a 6 DOF nonlinear aircraft model. 
The simulation results showed the potential of using this approach to achieve fault 
tolerant control. Based on this observation and some analysis, the paper proposes a 
multiple channel adaptive control scheme that makes use of the functionally redundant 
actuating and sensing capabilities in the model, and explains how to implement the 
scheme to tolerate actuator and sensor failures. The conditions, under which the 
scheme is applicable, are stated in the paper. 
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1. INTRODUCTION 

In applications where a system is expected to 
track a variety of reference trajectories, use of 
adaptive control becomes highly desirable. This 
is the case with highly maneuverable aircraft. In 
addition, one of the most thorny issues facing 
control engineers is the lack of an accurate model 
based on which a feedback compensation strategy 
necessary to guarantee the system performance is 
devised. The situation becomes more prominent 
when failures occur during a system operation 
such as that encountered by aircraft in a flight 
mission. Two approaches have been attempted 
to deal with the latter situation. One approach 
begins with a process of detecting and identifying 
parameter changes or faulty subsystems, which is 
then followed by a process of adjusting the control 
law. This approach sometimes suffers from high 
risk and low speed, attributed mainly to the de- 
tection and identification process. An alternative 
is to use adaptive control (Ahmed, et ah, 1991; 
Bodson, et ah, 1997; Chandler, et ah, 1995; & 
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Saeks, 1998). The reader is referred to Bodson 
(1997) for an overview of this approach to fault 
tolerant flight control. It is claimed in Saeks, Nei- 
dhoefer, Cox, and Rap (1998), with little analysis, 
that by using a multivariate variant of the Seraji 
(1989) decentralized adaptive control, originally 
developed for applications in the field of robotics, 
the following features can be attained: the design 
of a flight control law requires very little a pri- 
ori knowledge of the plant dynamics; the imple- 
mentation of the control law is straightforward; 
no explicit plant model identification is needed 
during adaptation; and the asymptotic stability 
of the flight control system is guaranteed. These 
all appear to be highly desirable for fault tolerance 
in flight control. 

In order to determine the ability of the decentral- 
ized adaptive control (Seraji, 1989) in handling 
initially poorly known aircraft models, especially 
aircraft models subject to changes due to unantic- 
ipated adverse operating conditions, and to assess 
the need for involving additional learning schemes 
such as neural networks, a study was conducted 
recently by Nikulin and Wu (1999), and Heimes 
(1999) using the NASA-Dryden 19-state nonlinear 
aircraft model (Brumbaugh, 1994) as the testbed. 



The model includes full envelope aerodynamics, 
propulsion system dynamics, actuator dynamics, 
and atmospheric model, for which a decentralized 
adaptive control law was attempted. The results 
of our study is reported briefly in Section 2. In 
Section 3, a multiple channel configuration is pro- 
posed for enhancing fault tolerance of the flight 
control system. 

2. DECENTRALIZED ADAPTIVE CONTROL 

There are five control inputs (aileron, symmet- 
ric stabilator, differential stabilator, rudder, and 
thrust) and nineteen measured states in the 
NASA-Dryden model (Brumbaugh, 1994). The 
model can be represented by a nonlinear state 
equation 

x = f(x , u) 

Our discussion on applying the Seraji control law 
will be focused on the linearized model of the form 

x = Ax + Bu, y =- Cx (1) 

at a particular operating condition. Since the Ser- 
aji control law is performance-based (as opposed 
to model-based), its dependence on the model 
fidelity is much less critical. The consideration on 
the effect of nonlinearity over the entire opera- 
tional envelope will be given at a later point. 

The Seraji control law, being fully decentralized, 
requires a one-to-one correspondence between the 
inputs and the outputs. Therefore, the dimension 
of y must be chosen to be equal to the dimension 
of u. Suppose dim(y)=dim(u) = m. The transfer 
function description of the linearized model is 
given by 

Y(s) = C[sl - A]~ l BU(*) = ^U(s) 

where d is a monic polynomial of degree and N 
is an m x m polynomial matrix. Consider s as a 
differential operator, which leads to the final form 
of what we call a design model 

d(s)y(t) = N(s)u(t) = T(t). (2) 

The design objective is to determine the control 
law u so that the asymptotic tracking of a selected 
set y of the measured states to a specified class 
of reference {y r e/} is achieved. The design model 
must be chosen carefully to satisfy the following 
properties. 

(a) y must have the same dimension as u. This is 
a constraint imposed by using the Saraji decen- 
tralized approach. 

(b) All interested, and potentially unstable states 


must be observable from y . This is a basic require- 
ment on fe edback variables of all control systems, 

(c) The transfer function from u to y must be min- 
imum phase, i.e., N~ l {s) must be a stable transfer 
matrix. This will guarantee that the transforma- 
tion from the ficticious control signal T to the 
actual control signal u is through a stable filter. 

Let r denote one of the components of T, p denote 
the corresponding component in y , and e = p re f~ 
p is the corresponding tracking error. Then we 
have the following. 

Theorem 1. The decentralized control law for (2) 
in one channel is given by 

v— 1 v 

T{t) = f(r) + Y + Y 0i(O/*Ji/(t).(3) 

1=0 1=0 

The auxiliary signal /, feedback gains and 

feedforward gains q^s are governed by the follow- 
ing equations 

/ = Sr + pr (4) 

k l ^a i (re^)+0 i (re {i) ) {l) ( 5 ) 

<1, = 7 i (r^e f)+ X i (V/e fY" ( 6 ) 

where < 5 , o n and 7 * are positive constants, p, i3i 
and A t are nonnegative constants, and 

r=p„ie-\ + /we ( " _1) . (7) 

p ul through p uv are from the last row of the 
solution P of the Lyapunov equation 

A t f P + PA f + Q = 0, Q = diag{Ki , •**, k„} > 0. (8) 

comes from the controller canonical realiza- 
tion Bf.Cf) of a desired v th order stable 
reference error filter, and positive numbers «i, K 2 , 

■ • - , k v must be chosen so that p vx > 0 V i, and 
Pvi > 0. The control law given by (3) through ( 6 ) 
ensures asymptotic reference tracking provided 
that p and its first u— 1 derivatives are slow signals 
relative to the auxiliary signal /. When this condi- 
tion is not satisfied, an additional constraint must 
be imposed in order to guarantee the asymptotic 
tracking. The constraint is given by 

t 

Mt)[vAt) ~ /( 0) -6 j r{a)da) < -^V'o(t), Vt, (9) 
0 

where 

Pd(t) = [d(s) - 

and Vo(t ) is the derivative of the Lyapunov func- 
tion used to derive the adaptation gains with Pd(t) 
set to zero, given by 

V V 

V 0 (t) = -->'Qe - 2r 2 {p + YMe (, '- l) )' 2 + Y X '^eff}- 

i— 1 i=0 



Vector e in the above expression is given by 

e = [e-e m (e - e m ) (1) ••• (e - e m ) (l/_1) ]', 

where e m is the state of the desired vth order 
stable reference error filter {Ap, Bp, Cp). □ 

The proof of this theorem draws heavily from the 
work of Seraji (1989). The modification of the 
Seraji control law to the current form (3)-(6) was 
suggested in Saeks (1998), where neither a proof 
nor conditions (7)--(9) were given. Since our proof 
for Theorem 1 is rather long, it is not included 
in this paper. The details of the proof will be 
provided upon request. 

Our design of the decentralized adpative control 
law for the NASA-Dryden model was successful. 
With the aid of a differential evolution strategy 
(Heimes, 1999) we are able to determine a single 
set of parameters in the adaptation laws (3)- (6) 
for all five flight conditions listed in Table 1. 
On the other hand, the design is by no means 
a straightforward matter. Many practical issues 
must be dealt with. The rest of the section is 
focused on discussing some of the issues: the 
model fidelity issue, the reference model issue, and 
the high order derivative issue. 


FI. cond. 

Altitude(ft) 

Mach 

Airspeed(ft/s) 

# 1 

17900 

0.62 

650 

# 2 

9800 

0.5 

540 

#3 

39800 

0.6 

580 

#4 

39800 

1.4 

1355 

# 5 

9800 

0.9 

973 


Table 1 


(A) Model fidelity. It can be seen that, besides the 
degree of d(s ), no other information about d(s) or 
N(s) is required to obtain control law (3). T(s), 
however, is only a fictitious control signal that 
cannot be applied to the aircraft. One must have 
the full information on iV, in order to calculate 
the real control signals 

u(s) = JV - 1 (s)r(s). (10) 

On the other hand, when the above input is 
applied 

«w = l r,wrw ’>' 

as if T were directly applied to plant I/d(s). 
Unfortunately, filter A r_1 (s) that comes from the 
design model generally does not match the N~ l (s) 
part of the true plant model. To distinguish a 
design model from a true model, subscript U D 
(design)” and subscript “T (true)” will be used 
wherever appropriate. One source of mismatch be- 
tween A r ^ 1 (s) and A'f 1 (s) is that the current sys- 
tem operating point may be different from the one 


at which the linearized design model is extracted. 
The seventy of mismatch varies from point to 
point within the operational envelope. Therefore, 
jV^ 1 ^), when cascaded to the input of the plant, 
introduces extra poles and cross-channel couplings 
because Nr(^)Np l ( 5 ) ^ /. For the NASA-Dryden 
model, A'^^s) is an extremely slow filter that 
it has a settling time of approximately 400 sec., 
which dominates the response of the system. To 
overcome this problem, we replaced low pass filter 
Np l (s) bj all pass filter A^ 1 (0). In doing so the 
residual signal [N^ l (s) - Np\0)]T(s) is ignored, 
in additional to the above mentioned mismatch. 
The residual signal is generally a fast signal con- 
taining the high frequency components of T, and 
thus easily violates the condition of Theorem 1. 
To satisfy die sufficient condition for Lyapunov as- 
ymptotical stability, a term — rqAy can be added to 
the right hand side of (5), called a rr-modification 
(Ioannou, 1986), as a more practical alternative 
to verifying inequality (9). As a result, a residual 
tracking error in the order of \fa is introduced. 
Cushioned by the sufficiency of the Lyapunov sta- 
bility theorem, the control system can sometimes 
be spared of instability without having to take 
this extra step, as we have observed in our design. 
However, our freedom in choosing parameters <5, p, 

■ • ■ in adaptation laws (4)-(6) becomes severely 
restricted. Suppose the extra step has been taken 
to guarantee the Lyapunov asymptotic stability, 
u(t) can be solved from 

N D (0)u(t) = T(t) 

The above 1 equation can be regarded as an effort 
to scale the ficticious control signals to the dimen- 
sions of the real input signals. Through the above 
equation c ross channel coupling is reintroduced. 

(B) Reference models. Much effort goes into the 
selection of the reference models (Heimes, 1999). 
The speed control reference model, lateral axis 
reference model, and longitudinal axis reference 
model are all selected are all selected to be suffi- 
ciently conservative to allow a single set of para- 
meters in the adaptive control law to work for all 
flight conditions, and yet sufficiently aggressive to 
allow a reasonable aircraft tracking performance. 

(C) High order derivatives. Because of the way the 
problem is formulated, high order derivatives with 
respect tc the tracking error and reference sig- 
nals enter the control law (3). Though derivative 
terms are desirable for improving the transient 
performance, by nature they are difficulty to be 
calculated accurately and they tend to amplify the 
noises. The extra room provided by the Lyapunov 
theory allows us to settle with a design that uses 
up to the third order derivative in any control 
channel, though the NASA-Drvden model could 




potentially lead to a 19th order derivative of the 
tracking error and of the reference signal in the 
control law. 


3. FAULT-TOLERANCE CONSIDERATION 

Our ultimate interest lies with the potential fault- 
tolerant capability of the decentralized adaptive 
control law. From the analysis in the previous sec- 
tion, the mild requirement on model fidelity of the 
Seraji control law motivates us to investigate the 
degree a control law can extend its authority in a 
secondary axis when loss of control effectiveness 
occurs in that axis. For example, the ailerons and 
the stabilators in the NASA-Dryden model can 
be made to provide redundant control authorities 
in the roll, and the rudder and the engines can 
be made to provide redundant control authorities 
in the yaw axis. System instability, reduced ma- 
neuverability, and cross-coupling are among the 
consequences of sensor /actuator failures and con- 
trol surface damage (Bodson, 1997). Therefore, 
it is important that a multivariate approach be 
taken to take advantage of the existing functional 
redundancy. 

Fig.l shows the measured roll/yaw rates and the 
control signals of a simulated seenario-a coordi- 
nated turn, carried out at an altitude of 17900 ft 
and an airspeed of 650 ft/sec using the stabilators 
and the rudder of the NASA-Dryden model. The 
decentralized adaptive control described in the 
previous section is used. Fig. 2 shows the same set 
of signals resulting from the same maneuver at the 
same flight condition. This time however, an 80% 
loss of right stabilator effectiveness is simulated. 
In this case, only the auxiliary control signal (4) is 
applied, which represents the conventional portion 
of the adaptive control law (4)-(6). It can be seen 
that the aircraft performance is severely degraded. 
After turning on the adaptive control law, as 
shown in Fig. 3, the system is able to maintain the 
trim and restore the tracking under the failure 
condition. It is also observed (but not shown in 
the plots) that the side slip angle in this situation 
is kept within the prescribed limit of ±2°. The 
auxiliary signal, and the adaptation gains for the 
roll rate tracking error and the roll rate reference 
signal are shown in Fig. 4. This indicates that with 
adequate redundancy in the plant it is possible 
to design an adaptive control law that tolerates 
severe adverse conditions. The rest of the paper 
devotes to exploring such potentiality through a 
new fault tolerant configuration that requires a 
minimum local monitoring effort at the sensors 
and the actuators. Investigation is on going to 
implement the new configuration through simu- 
lations. Fig.l through Fig.4 are all shown in the 
last page of the paper. 


Suppose a design model at a particular operating 
point takes the form 

i* = Ax f B\U\ T 13 2 U‘2 j y\ — C\x* U2 — C2*r(H) 

where .r 6 R n , u x G R m] , u 2 € y x € FT\ 

and ?/2 € R rn ' 2 . Typically in this kind of applica- 
tions mj m 2 — 3. The division of the inputs is 
based on the result of an analysis so that each n t 
can provide adequate full control authority, which 
requires sufficient controllability of (A, B x ). y\ and 
y 2 are two alternatives, each adequately reflecting 
the requirements on the tracking performance, 
which requires sufficient observability of (Cj.d). 
A combined measure of controllability and observ- 
ability, called by Wu, Zhou, and Salomon (2000) 
a reconfigurability, should be imposed. More than 
two sets of inputs and outputs should be consid- 
ered whenever the redundancy allows. In addition, 
the input to output transfer matrix defined by 
(variable 7 is suppressed) 


Nil A 12 


7/1 


A r 2i A 22 


u x 

N 

771 

_ 7/2 , 


d 

u 2 

“ 7 



is a minin um phase system, i.e., N~ l (s) is a sta- 
ble transfer matrix. We also require that N~ 1 (s)’s 
are stable transfer matrices. d(s) is a monic poly- 
nomial of degree v in s, N(s), and A r ^ (s)’s are 
mi 4- 7712 x mi +77i2, and rn i x m, 3 polynomial ma- 
trices, respectively, of degrees less than v. Again 
consider s as a differential operator, and define 
the design model 

d(s)yi{t) — jVn(.s)wi(t) + iV 12 (s)u 2 («) = 7\(/)j(13) 
d(s)y 2 (t) = 1 V 21 («)«,(*) + N 22 (s)u 2 (t) = T 2 (t)(l 4) 

We are now in a position to state the asymptotic 
tracking problem. Determine T\ and T 2 so that 

2/1 0 ) ( “+° V\ref, andy 2 (0 f -7° 2/2re/- (15) 

Obviously, if T\ and T 2 can be produced satisfac- 
torily, u\ and u 2 can be obtained by sending T\ 
and T 2 through the stable filter A r_1 (s). Moreover, 
the individual components of yi and those of T x 
have one t o one correspondence. 

Theorem 2. Suppose T x and T 2 in (13) and 
(14) solve the asymptotic tracking problem. When 
some of the components of y x become unavailable 
(= 0), control law Tj , j ^ i solves the asymptotic 
tracking problem. When some of the components 
of U{ become unavailable (=0), control law u 3 — 
N ~ 1 Tj , j / j solves the asymptotic tracking 
problem. 

Proof. Suppose some components of yt are no 
longer available. From (3) it is seen that Tj , j ^ i 
depends only on y 3 and y Jre /. Therefore y 3 — ► 



ijjref , provided that (9) is satisfied for design 
model d(s)yj{t) = Tj(t). Since v/ 7 has been se- 
lected to reflect the tracking requirement, T 3 
solves the asymptotic tracking problem. Suppose 
now that u t becomes unavailable. 

Tj = A ij u x + NjjUj =- A jjUj. 

Again Tj, j ^ i solves the asymptotic tracking 
problem, i.e., y 3 — ► i/j re /, provided that (9) is 
satisfied for design model d{s)y 3 ,t) — Tj(£). Since 
Uj (s) = A T “ 1 (s)T J (s), and A'" 1 is stable, u 3 solves 
the asymptotic tracking problem. □ 

Note that if model fidelity is an issue, ^-modification 
(Ioannou, 1986) may be required in order to 
achieve asymptotic tracking. Fig. 8 shows one im- 
plementation among many possible implementa- 
tions of the multiple channel configuration. This 
implementation requires the local monitoring of 
sensors and actuators, w r hich is a much less de- 
manding task than the task of diagnosis. In this 
setup, when measurement y\, or input u\ be- 
comes faulty, filter A^ 1 is set to zero. When 
measurement y i, or input u\ becomes faulty, filter 
(A 22 - A r 2 i A^'V A r i 2 ) -1 is set to zero. In general 
however, no reconfiguration is required as long 
as the remaining control input set contains at 
least all components of u t , or that the remaining 
measurement set contains at least all components 
of in addition to satisfying the conditions of 
Theorem 1. 



An Implemen tation of Filter N ’ 1 


(b) 

Fig. 5 A multi-channel configuration 

Normally a better tracking performance can be 
expected if the parameters in the adaptation law 
(4)-(6) are scheduled for different flight conditions 
(Nikulin and Wu, 1999). When information is 
available, the parameters of filter N _1 can also be 
scheduled. In that case, the single filter in Fig.5 is 
replaced by a bank of filters scheduled for various 
flight conditions. 

4. CONCLUSIONS 

This paper presents some results of our initial 
investigation on the feasibility to extend a decen- 


tralized adaptive control scheme (Seraji, 1989) to 
flight control. Such an extension is suggested in 
Saeks (1998). The paper examines the conditions 
under which the decentralized adaptive control 
can be applied, and explicitly implements such a 
decentralized adaptive control law on the NASA- 
Drvden model (Brumbaugh, 1994). Though the 
analysis presented in the paper validates the de- 
centralized control law, there is no easily checkable 
conditions for guaranteeing the asymptotic track- 
ing. Therefore, extensive simulation plays a very 
important role in ensuring a successful design. Our 
simulation and analysis results have shown the 
promise ol the method in handling severe aircraft 
impairment, and severe under modeling. Based 
on this observation, a multiple channel configura- 
tion is proposed for achieving fault-tolerance. This 
configuration requires the implementation of a full 
scale multivariable design which is being pursued. 
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Fig.l Measured roll/yaw rates and Fig. 3 Measured roll/yaw rates and 

rudder/stabilator control signals rudder /stabilator control signals 


80% Loss of Stabliator Effectiveness 
(conventional control) 
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Fig. 2 Measured roll/yaw rates and 
rudder/stabilator control signals 


80% Loss of Stabliator Effectiveness 
(decentralized adaptive control) 
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Fig.4 Stabilator auxiliary signal, and 
leading feedback, and forward adaptation gains 










